First-principles predictions of low-energy phases of multiferroic BiFeC>3 
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We used first-principles methods to perform a systematic search for potentially-stable phases of 
multiferroic BiFeC>3. We considered a simulation cell compatible with the atomic distortions that 
are most common among perovskite oxides, and found a large number of local minima of the energy 
within 100 meV/f.u. of the ferroelectric ground state. We discuss the variety of low-symmetry 
structures discovered, as well as the implications of these findings as regards current experimental 
(e.g., on thin films displaying super-tetragonal phases) and theoretical (on models for BiFeOa's 
structural phase transitions) work on this compound. 
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I. INTRODUCTION 

Perovskite oxide BiFeC>3 (BFO) continues to reveal it- 
self as one of the most intriguing materials of the day. 
Not only does it remain the most promising magnetoelec- 
tric multiferroic for applications at room temperature, 
but it also has been shown recently to display a variety 
of novel fundamental effects.^ Such findings range from 
an increased conductivity at specific ferroelectric domain 
walls^ to new structural phases in thin films^^ with po- 
tentially useful response properties!^ 

The present work originated from our on-going 
research on enhancing the properties of BFO by 
forming solid solutions such as BiFei-^Co^OgP and 
Bii-zRE^FeOa with RE = La, Sm, DyP While inves- 
tigating the chemically induced structural transitions, it 
became clear we needed to have a thorough and unbiased 
strategy to search for possible structural phases beyond 
those reported in the literature. Interestingly, when we 
applied such a scheme to BFO itself, we found plenty of 
low-symmetry phases that are local minima of the en- 
ergy. Here we describe the lowest-energy structures that 
we discovered, i.e., those most likely to be observed ex- 
perimentally. We discuss the origin of the large variety 
of distortions found in the calculations, and the possibil- 
ity of capturing BFO's structural richness within simple 
models. Further, we comment on the implications of our 
findings as regards current experimental work on BFO in 
both bulk and thin film forms. 



II. METHODOLOGY 

For the simulations we used the local density (LDA^ 
and generalized gradient (PBB^S! and PBEsoP^) approx- 
imations to density functional theory (DFT) as imple- 
mented in the VASP package. 12 ' A "Hubbard- U' scheme 
with U = 4 eV was used for a better treatment of iron's 
3c? electrons,^ the corrected DFT functionals will thus 
be referred to as LDA+ U, PBE+U, and PBEsol+U. We 
used the "projector augmented wave" method to repre- 
sent the ionic coresP solving for the following electrons: 
Fe's 3s, 3p, 3d, and 4s; Bi's bd, 6s, and 6p; and O's 




(b) (c) 

FIG. 1. (Color online.) (a) 40-atom supercell of BiFeOs (ex- 
tra periodic images of some O and Bi atoms have been in- 
cluded for easier visualization). The O atoms occupy the 
vertices of the octahedra plotted, which contain Fe atoms at 
their centers; the rest of the atoms are Bi. (b) Same cell 
as (a), illustrating a C-AFM arrangement of Fe atom spins, 
(c) Same as (b), but with a G-AFM spin arrangement. All 
the phases considered in this work have unit cells that are 
distortions of the one depicted here. 

2s and 2p. (We checked that qualitatively correct re- 
sults can be obtained without considering semi-core elec- 
trons.) Wave functions were represented in a plane- wave 
basis truncated at 500 eV, and a 2x2x2 fc-point grid was 
used for integrations within the Brillouin zone (BZ) cor- 
responding to the 40-atom cell of Fig. [T] The calculation 
conditions were checked to render converged results. 
We worked with the 40-atom cell depicted in Fig. [T] 
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which is obtained by doubling the 5-atom cell of the ideal 
perovskite structure along the three Cartesian directions, 
denoted by x, y, and z in the following. This cell is 
compatible with the structural distortions that charac- 
terize the low-symmetry phases of many perovskite ox- 
idesP^ (1) ferroelectric (FE) patterns associated with ir- 
reducible representation (symmetry labels correspond 
to the BZ of the 5-atom cubic cell); (2) anti- ferroelectric 
(AFE) modes associated with zone-boundary q points 
(Jf-like, M-like, and R); and (3) anti-ferrodistortive 
(AFD) patterns corresponding to any combination of in- 
phase (M^) and anti-phase (R±) rotations of the C>6 
octahedra around the Cartesian axes. This cell is also 
compatible with the anti-ferromagnetic (AFM) spin ar- 
rangements known to be most relevant for BFO, i.e., the 
C-AFM and G-AFM orders sketched in Figs, [ljb) and 
[ljc), respectively. 

To explore all these possibilities we considered a large 
number of starting configurations for our structural relax- 
ations. Specifically, we considered: (1) all AFD patterns 
consisting of either an in-phase or an anti-phase rota- 
tion around each Cartesian axis (i.e., those expressible in 
Glazer's notation^; (2) various FE and AFE patterns 
constructed by off-centering the Bi cations; (3) cells with 
cubic, tetragonal, and orthorhombic shapes; (4) G- and 
C-AFM orders as well as a few attempts with other spin 
arrangements. This added up to more than 300 starting 
configurations. In all cases, we first ran a short molec- 
ular dynamics simulation with random initial velocities 
(thus breaking all symmetries) , and then performed a full 
structural relaxation. We used the PBE+ U functional for 
this structural search. The lowest-energy configurations 
obtained were confirmed to be minima by checking their 
stability against ionic and cell distortions. 



III. RESULTS 

A. Lowest-energy phases found 

Our search led to a wealth of local minima with ener- 
gies in a range up to 200 meV/f.u. above BFO's ground 
state. Table I lists the lowest-lying solutions; we show 
their PBE+U energy directly obtained from our struc- 
ture search, as well as the energies obtained by relaxing 
the PBE+ U structure using the PBEsol+ U and LDA+ U 
functionals. Note that the energy differences between 
phases are strongly dependent on the DFT functional; we 
will address this issue below. Table I also includes a short 
description of the phases found, which we label by their 
atomic space group and type of AFM order (e.g., i?3c-G 
for the ground state); the complete structural informa- 
tion and computed polarization values^ are given in Ta- 
bles II and III. Let us note that our work with BFO and 
other compounds confirms that PBEsol is more accurate 
than PBE and LDA for predicting the atomic structure of 
individual phases.^ Thus, the crystallographic data re- 
ported here correspond to PBEsol+ U-relaxed structures. 



Finally, Fig. 2 shows sketches of the structures obtained, 
and the most relevant distortion modes are depicted in 
Fig. 3. 

All the functionals correctly predict the i?3c phase with 
G-AFM spin order as the ground state of BFO. This 
phase displays a spontaneous polarization along the [111] 
Cartesian direction, and anti-phase Oq rotations around 
the same axis (a~a~a~ in Glazer's notation). 

We also found two orthorhombic phases that are sim- 
ilar to i?3c-G in that they involve a relatively small dis- 
tortion of the ideal cubic cell and favor the G-AFM order: 
Pnma-G and Pna2\-G. 

The Pnma-G structure is paraelectric (PE). As shown 
in Table I, it is characterized by an Oq rotation pat- 
tern (a~a~b + ) that involves anti-phase rotations around 
[110] and in-phase rotations around [001]. This phase is 
the ground state of many perovskites, LaFeOa being the 
most relevant one for the current discussion. Interest- 
ingly, BFO's Pnma-G phase can be aptly described as 
AFE, because the Bi cations present large anti-polar dis- 
placements in the (001) plane (associated with the 
mode of Table I and Fig. 3) ; the computed off-centering 
of the Bi cations is about 0.3 A. (Such an AFE pattern 
is allowed by symmetry in LaFeOa too; in that case we 
obtain La off-centers by about 0.2 A.) 

The Pna2i-G phase is similar to Pnma-G, but with 
an additional FE distortion along the axis of the in-phase 
rotations. As compared with that of Pnma-G, the 40- 
atom cell of the Pna2\-G phase is elongated along the 
polarization direction; this reflects the usual coupling be- 
tween the FE distortion and strain observed in perovskite 
oxides. 

Regarding magnetism, the R3c-G, Pnma-G, and 
Pna2i-G phases display strong AFM exchange couplings 
between neighboring Fe ions, as evidenced by a large en- 
ergy splitting of more than 200 meV/f.u. between the 
G-AFM and ferromagnetic (FM) configurations. This is 
consistent with the high magnetic ordering temperature 
observed in bulk BFO. 

In addition, we found a number of phases that involve 
a large stretching of the ideal cubic cell along the z di- 
rection, with c/a aspect ratios approaching 1.3. In the 
following we will generically refer to such structures as 
super-tetragonal or T phases. They all favor the C-AFM 
order (see Fig. 1), the parallel spin alignment occurring 
along the stretched lattice vector. The magnetic interac- 
tions along z are weak, as evidenced by an energy split- 
ting of about 5 meV/f.u. between the C- and G-AFM 
orders; accordingly, the ordering temperatures will be rel- 
atively low. Three of these phases are monoclinic (Cc-C, 
Cm-G, and Pc-G) and one is orthorhombic (Pna2i-G); 
all of them are ferroelectric with a very large polarization 
component along [001] (see computed values in Table II). 

More specifically, the Cc-C phase presents a polariza- 
tion in the (110) plane, as well as relatively small AFD 
distortions. This type of monoclinic phase is usually 
termed a similar phase has been studied theoret- 

ically in connection with the super-tetragonal structures 
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TABLE I. Energies and distortions of the most stable energy minima found, as well as a few saddle points (six bottom phases) 
included for reference. Columns 2-4: Energies obtained with different DFT functionals. Note Pna2\-G goes to Pnma-G 
when relaxed with PBEsol+U and LDA+U. Columns 5-8: Distortions from the ideal cubic perovskite structure (Pm3m) that 
characterize the phases. In all cases the FE and AFD modes fully determine the symmetry breaking. A generic [x, y, z] 
FE (AFD) distortion involves displacements (06 rotations) along (around) the x, y, and z Cartesian axes. We indicate the 
dominant FE and AFD distortions in bold. Column 8 includes other modes with a significant contribution (at least 10% of 
largest one). The mode analysis was done with the isodisplace softwarep' note that q-points indicated in symmetry labels 
constitute default choices and do not always correspond to the actual distortion modulation (e.g., the X$ and X^ AFE modes 
in the Table are actually modulated along the z direction). 
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TABLE II. Computed PBEsol+U lattice parameters (corresponding to the 40-atom cell of Fig. 1) and polarization values for 
the six stable phases of BFO listed in Table I. The polarization direction is given in a Cartesian reference that corresponds 
almost exactly with the 40-atom cell vectors. For comparison, we also include the result for the P4mm-C structure. 



Lattice parameters Polarization 



Phase 


a (A) 


b(k) 


c(A) 




PC) 


7(°) 


Magnitude (C/m 2 ) 


Direction 


Pc-C 


7.500 


7.500 


9.489 


88.1 


88.1 


89.7 


1.20 


(0.29, 0.29, 0.92) 


Cm-C 


7.380 


7.608 


9.533 


86.6 


90.0 


90.0 


1.50 


(0.00, 0.30, 0.95) 


Pna2i-C 


7.515 


7.515 


9.452 


90.0 


90.0 


90.0 


1.39 


(0.00, 0.00, 1.00) 


C*c-C 


7.527 


7.527 


9.444 


88.0 


88.0 


90.0 


1.45 


(0.23, 0.23, 0.94) 


Pnma-G 


7.830 


7.830 


7.770 


90.0 


90.0 


87.6 







P3c-G 


7.893 


7.893 


7.893 


89.5 


89.5 


89.5 


0.91 


(0.58, 0.58, 0.58) 


P4mm-C 


7.414 


7.414 


9.526 


90.0 


90.0 


90.0 


1.52 


(0.00, 0.00, 1.00) 



observed experimentally in BFO mms p512LLSK' 

The Pc-C phase is very similar to Cc-C as regards the 
polar distortion (thus, it is also Ma), but it displays a 
different Og-rotation pattern. 

The Cm-C phase displays a polarization in the (100) 
plane, and the cell is significantly distorted in the xy 
plane; such a monoclinic phase is termed Mc- 

The Pna2i~C phase is very similar to the Pna2 1 -G 
structure discussed above, the stretching of the cell and 
development of polarization coinciding with the axis of 
the in-phase rotations. 

Note that all these T phases can be viewed as distorted 
versions of the ideal super-tetragonal P4mm-C structure 



listed in Table I. Interestingly, we found that this PAmm- 
C phase, which is the ground state of BiCo0 3 plis a saddle 
point in BFO's energy landscape. 

Our results thus reveal an intricate energy landscape, 
especially regarding structures with a large c/a ratio. 
In this sense, it is interesting to note that some of the 
phases reported here are small distortions of higher- 
symmetry structures. For example, the Cm-C phase can 
be shown to be a Pm-C structure distorted by the 
[0, y, 0] mode listed in Table I; by moving from the Pm-C 
saddle point to the Cm-C minimum, BFO gains about 
1 meV/f.u. Similarly, the reported Pc-C phase is con- 
nected with a higher-symmetry Cm-C structure via a 







ft T • ° 



(a) Pc-C 



(b) Cm-C 






(c) Pna2 1 -C 



(d) Cc-C 





(e) Pnma-G 



(f) R3c-G 



FIG. 2. (Color online.) Energy minimum configurations obtained, (a)-(d) C-AFM super-tetragonal phases; in the left (right) 
image the c axis is perpendicular (parallel) to the page, (e)-(f) G-AFM phases; two pseudo-cubic axis are equivalent in (e), 
with the left (right) figure having the non-equivalent axis perpendicular (parallel) to the page; the three pseudo-cubic axis are 
equivalent in (f). The atomic species can be identified as in Fig. [I] 



M^-[0, 0, z] distortionPsl Given BFO's manifest complex- 
ity, we tend to view the phases of Table I as a probably- 
incomplete list, just indicative of the rich variety of stable 
structures that this compound can present. 

Finally, let us stress we explicitly checked that all the 
above phases are local minima of the energy, a fact that 
is remarkable since some of them (e.g., the pairs formed 
by Pnma-G and Pna2i-G, or Cc-C and Pc-C) are rather 
close structurally. It is also interesting to note that mon- 
oclinic phases with such small primitive cells may be en- 
ergy minima by themselves, i.e., without the need of any 
stabilizing (electric, stress) fields. Note that, to the best 
of our knowledge, monoclinic phases in bulk perovskite 
oxides tend to be associated with complex solid solutions 
or large unit cells. Examples of the former are the mon- 
oclinic Ma phase that occurs in prototype piezoelectric 



PbZri-^TizOaP*! and the monoclinic Mc phase of rc- 
laxor PbZn 1/3 Nb 2/ 303-PbTi0 3 P Examples of the lat- 
ter occur in BiMnC>3 and BiScOa; see the discussion in 
Ref. [STJ It was thus unexpected to discover that bulk 
BFO presents such a collection of simple low-symmetry 
minima of the energy. 



B. Energy differences between phases 

The relative stability of the phases discussed above is 
quantified by the energy differences between them. Dis- 
turbingly, Table I shows that such energy differences are 
strongly dependent on the DFT functional used to com- 
pute them. By switching functional we obtained changes 
in relative stability - e.g., Pnma-G is more stable than 
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FIG. 3. Illustration of atomic displacements for different sym- 
metry modes of BFO: (a) soft FE mode, (b)-(f) secondary 
modes mentioned in Table I. Only displacement directions, 
not magnitudes, are indicated; for the (a) case, the PBEsol+fJ 
computed atomic displacements are quoted in the caption of 
Fig. [7] White, grey, and black circles represent Bi, Fe, and O 
atoms, respectively. 



the T phases according to PBEsol+ U and LDA+ U, but 
less stable according to PBE+ U - and even the loss of 
stability of one phase - i.e., the Pna2i-G phase is sta- 
ble for PBE+ 17, but the relaxation of this structure with 
PBEsol+I7 and LDA+[7 leads to the Pnma-G solution. 
Thus, we have to ask ourselves: Does any of these func- 
tionals provide an accurate picture of the relative sta- 
bility of BFO's phases? Noting that PBEsol is gener- 
ally more accurate regarding the structural description 
of individual phases, can we just rely on the PBEsol+tJ 
results? 

One would like to address this issue by resorting to a 
higher-level first-principles theory. However, performing 
quantum Monte Carlo calculations, which are the refer- 
ence for accuracy in this context, is well beyond the scope 
of this work. Simpler schemes like the so-called hybrid 



TABLE III. Energy minima structures of Table I as obtained 
from PBEsol+U calculations. In the case of the Pna2\-G 
phase, the PBE+ U result is given (see text). 

Pc-C a = 7.291 A b = 5.291 A c = 5.315 A 
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Bi 


2a 
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0.9418 


Fe 
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c = 9.452 A 




a 


= p = 7 = 90° 






Atom 


Wye. 


X 


y 


z 


Bi 


4a 


0.5451 


0.4799 


0.4590 


Fe 
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Cc-C 
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0.4829 
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Fe 
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0.2630 
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0.3290 
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4a 


0.3405 
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0.4593 



functionals, which are usually considered to be more ac- 
curate than DFT for insulators like BFO, are not well 
tested for quantifying relative stabilities in cases like this 
one. Moreover, structural predictions with hybrids have 
been shown to depend strongly on the underlying gener- 
alized gradient approximation,^ which invalidates them 
for the present purposes. 

Nevertheless, we were able to make a couple of mean- 
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TABLE III. (contd.) 



Pnma-G 


a 


= 5 650 A b — 


7.770 A 


c = 5.421 A 






a 


= P = 7 = 90° 






Atom 




Wye. 


X 


V 


z 


Bi 




4c 


0.0523 


1/4 


0.0100 


Fc 




4h 








1/2 







4c 


0.9722 


1/4 


0.5946 







8d 


0.2998 


0.0461 


0.3037 


R3c-G 






a = b = 


5.559 A c 


= 13.782 A 








a = = 


90° 7 = 


120° 


Atom 




Wye. 


X 


y 


z 


Bi 




6 a. 








0.0000 


Fc 




6 a 








0.7236 







18b 


0.3156 


0.2294 


0.1238 


Pna2\ 


G 


CI 


= 5.702 A b = 


5.507 A 


c = 8.036 A 


(PBE4 


-U) 


C\ 


— B — ■y — 90° 






Atom 




Wye. 


X 


V 


z 


Bi 




4a 


0.4435 


0.0016 


0.2194 


Fc 




4a 


0.5015 


0.5007 


0.4943 







4a 


0.2137 


0.7074 


0.0519 







4a 


0.1848 


0.6876 
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ingful comparisons with experiment. First, we studied 
the transition between the R3c-G and Pnma-G phases 
that is known to occur under hydrostatic pressure!^ 
We obtained (see Fig. 4) transition pressures of about 
2 GPa for LDA+U, 3 GPa for PBEsol+U, and 5 GPa for 
PBE4- U. Room-temperature experiments by Haumont et 
alW^ showed that at 3.5 GPa the i?3c-G phase trans- 
forms into a monoclinic C2/m structure with a large cell 
(made of 12 formula units), and that a second transi- 
tion at 10 GPa leads to the Pnma-G phase. These re- 
sults suggest that the i?3c-G and Pnma-G phases revert 
their relative stability at a pressure between 3.5 GPa and 
10 GPa, a bracket that can be shifted to 5-14 GPa if 
the transition lines are extrapolated to KP Thus, this 
comparison seems to indicate that the PBE+U is the 
most accurate theory for relative stability calculations, 
and that the LDA+ U should not be used for these pur- 
poses. We have reached similar conclusions in our work 
with Bii_a;La s Fe03 solid solutions^ in that case, the 
LDA+ U predicts a R3c-to-Pnma transition for a La con- 
tent that is clearly too small to be compatible with the 
experiments. 

Second, we computed the relative stabilities of these 
phases as a function of an epitaxial strain correspond- 
ing to a square substrate in the (001) plane, so as to 
determine the lattice mismatch needed to stabilize the 



PBE+U PBEsol+U LDA+U 




FIG. 4. (Color online.) Energy versus volume curves for the 
most stable phases of BFO. The labels at the top indicate the 
DFT functional used. The transition pressures mentioned in 
the text were obtained by computing the slope of the common 
tangent of the R3c-G and Pnma-G curves. 

large-c/a structures.^ As shown in Fig. 5, we obtained 
strain values of -2.3%, -4.0%, and -4.5% for PBE+U, 
PBEsol+U, and LDA+U, respectively. Experimentally 
it is known that a BFO-(001) thin film grown on SrTiOa 
(—1.5% misfit strain) displays a monoclinic structure 
that is an epitaxially-distorted version of the i?3c phase 
(such a phase is believed to be monoclinic Ma with the 
Cc space groupP^); we will denote this phase by R in the 
following. In contrast, when LaA103 substrates (—4.8% 
misfit strain) are used, a super-tetragonal T phase whose 
symmetry remains unclear,^ or a co-existence of the R 
and T phases^ has been observed. These results suggest 
that the energies of the R and T phases cross at an epi- 
taxial compression close to —4.8%. Hence, according to 
this criterion, and assuming that our large-c/a phases are 
good candidates to be the observed T phase, the PBE+U 
curves would be the least reliable ones. We have reached 
similar conclusions in our work with BiFei-^Co^Oa solid 
solut ions pSl where PBE+ U predicts an R-to-T transition 
for a Co content that is too small to be compatible with 
experiment. Further, these observations seem consistent 
with a well-known failure of the PBE approximation: it 
tends to render too large tetragonal distortions in ferro- 
electric perovskites! 31 * 37 * 

In conclusion, while the PBE+U and LDA+U ap- 
proaches seem to be rather accurate in some cases, they 
also render clearly wrong predictions in others. In this re- 
spect, PBEsol+ U seems to be a reasonable compromise, 
as it constitutes the overall most accurate DFT theory 
available to us. Nevertheless, because PBE+ U performs 
well as regards the relative stability of the i?3c-G and 
Pnma-G phases, we believe that the PBE+U prediction 
of the new ferroelectric phase Pna2\-G, structurally very 
similar to Pnma-G, deserves some attention. Finally, let 
us note that the choice of U also has en effect on the 
energy differences of Table I. Yet, for U values in the 3- 
5 eV range, such effects are small as compared with the 
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FIG. 5. (Color online.) Energy of various BFO phases as 
a function of the misfit (epitaxial) strain corresponding to a 
square (OOl)-oriented substrate. The labels at the top indicate 
the DFT functional used. Note that the R3c-G phase reduces 
its symmetry to Cc-G in these epitaxial conditions. 



ones we have discussed. 



IV. DISCUSSION 

Our results have direct implications for current ex- 
perimental work on the structural characterization and 
phase transitions of BFO, especially regarding the epitax- 
ially compressed films in which super-tetragonal phases 
were discovered. Further, they also provide us with in- 
formation that is relevant to the effective modeling of 
BFO's structural transitions, at both the macroscopic 
(Landau-type theories) and atomistic (effective Hamilto- 
nians) levels. In the following we discuss all these aspects. 
To conclude this Section, we comment on Bi's ability 
to form very different and stable coordination complexes 
with oxygen, as this seems to be the factor responsible 
for the observed structural richness of BFO. 



A. Implications for experimental work 



Our present results show that there are many possible 
T phases - e.g., the low-energy Pc-C, Cm-C, Pna2 1 -C 1 
and Cc-C structures that we found - that might cor- 
respond to the one experimentally realized in the BFO 
films. Indeed, as shown in Fig. 5, all our large-c/a phases 
are essentially degenerate in energy for values of the epi- 
taxial strain corresponding to a LaA103 substrate. More- 
over, at the PBEsol+ U level - which we have adopted as 
the DFT flavor of choice for BFO -, all these phases have 
their energy minimum at a misfit strain of about —4.8%, 
implying that any of them can form a stable BFO film 
under such epitaxial conditions. 

Because our T phases are an almost perfect epitaxial 
match with the LaA103 substrate, the structural and po- 
larization data in Tables II and III can be compared with 
the experimental results directly. Most remarkably, our 
results show that phases with very similar c/a ratios can 
display rather different polarization values. Indeed, the 
Pc-C phase (with a c/a of 1.27) presents P z as 1.1 C/m 2 , 
while the Cm-C and Pna2\-C phases (with c/o's of 1.26 
and 1.25, respectively) present P z as 1.4 C/m 2 . Hence, 
our Pc-C structure seems to be the best candidate to rep- 
resent the T phase realized in the BFO films investigated 
experimentally; the quantitative disagreement between 
the measured and predicted P z 's would be below 40%, 
a clear improvement upon the previously reported 90% 
difference. 

Let us also note that, because our T phases are so 
close in energy, the question of which one is realized 
experimentally may depend on subtle details not con- 
sidered in this work. Thus, for example, two of these 
phases (Pc-C and Cm-C) present no tilts (i.e., rotations 
around the [100] and [010] axes) of the 06 octahedra, 
which may make them preferable if the BFO films are 
grown on (001) substrates that clamp such distortions 
strongly. Similarly, a rectangular substrate might favor 
the Cm-C phase, whose cell tends to distort in the xy 
plane, etc. 

Fin ally, w e have very recently become aware of new 
resultaESEE] showing that both Mq and Ma monoclinic 
phases with large-c/a ratios can be realized in epitaxi- 
ally compressed BFO-(OOl) films. Such findings further 
support the physical relevance of the present study. 



1. Super-tetragonal phases in BiFeOz films 

The recent works by Bea efc aiPand Zeches et aiPhave 
shown that it is possible to obtain a novel phase of BFO 
if thin films are grown on strongly compressive substrates 
like LaAlO3-(001). Experimentally, this T phase presents 
a very large c/a ratio of about 1.23, and an out-of- plane 
polarization P z as 0.8 C/m 2 . First-principles studies^^D 
have identified the T phase with a monoclinic Cc struc- 
ture for which LDA+ U calculations predict c/a as 1.23 
and P z as 1.5 C/m 2 . Thus, there is a large quantitative 
discrepancy between theory and experiment as regards 
the value of P z , which suggests that the identification of 
the simulated and experimental phases may be incorrect. 



2. Structural transitions in bulk BiFeOz 

Our calculations were restricted to the limit of low 
temperatures, and do not allow for a conclusive dis- 
cussio n of temperature-driven effects and transitions in 
BFO 1221221 Nevertheless, a few comments can be made 
based on the obtained (large) energy differences between 
some relevant phases. Inde ed, our results seem consis- 
tent with experiment j 1 ! 24 ! 2 ^ showing that, as a function 
of increasing temperature, BFO's ferroelectric P3c phase 
transforms into an orthorhombic Pmna structure at 
T as 1100 K, to then become cubic Pmlm at T as 1200 K. 
More specifically, the PBEsol+ U results of Table I show 



that the i?3c-G and Pnma-G minima are very close in 
energy and constitute strong instabilities of the proto- 
type Pm3m structure, which lies about 900 meV/f.u. 
above them, as consistent with the fact that BFO's cubic 
phase can be observed only at very high temperatures. 
Moreover, R3c-G and Pnma-G constitute BFO's most 
stable phases, with a large margin over other structures 
(e.g., the ferroelectric R3m-G and Amm2-G, or para- 
electric R3c-G and IA/mcm-G, listed in Table I) that 
are common among perovskite oxides. Hence, our results 
seem incompatible with the R3c — > 1 4/ mem — > PmZm 
transition sequence obtained by Kornev et aiPH from 
Monte Carlo simulations of first-principles-derived effec- 
tive Hamiltonians; we found that the I A/ mem structure 
has a relatively high energy and is thus unlikely to occur 
instead of Pnma. 

As regards pressure-driven transitions, our results con- 
firm that under compression BFO's i?3c-G phase loses 
stability in favor of the Pnma-G structure! 27 * 38 ! Addi- 
tionally, it is worth noting that, at the PBE+J7 level, 
we found a Pna2\-G phase (see Table I) whose stabil- 
ity is also favored by compression and which nearly be- 
comes the ground state in the pressure range in which 
i?3c-G and Pnma-G revert their relative stability (re- 
sults not shown here). Given that PBE+I7 seems the 
most accurate DFT flavor for the description of these 
pressure-induced transformations (see Section III.B), it 
seems wise to bear in mind the possibility that such a 
Pna2\-G structure might occur, especially considering 
that the nature of the phase intermediate between RZc- 
G and Pnma-G remains unclear 



B. Implications for modeling work 

Our results clearly demonstrate that, in spite of its 
apparent simplicity, BiFe03 is extraordinarily complex 
from the structural point of view. In the following sec- 
tions we will quantify such a complexity, adopting the 
perspective of someone who is interested in determin- 
ing the simplest possible model, either macroscopic or 
atomistic, that captures accurately BFO's structural di- 
versity. Our analysis shows that BFO is much more chal- 
lenging to model than traditional ferroelectric perovskites 
like BaTi0 3 , PbTi0 3 , or even PbZr^Ti^Oa. 



1. Primary and secondary distortions in BiFeOs, 

By analyzing the BFO phases described in Table I, 
it is possible to identify three primary distortion types 
(or primary order parameters) whose occurrence can ex- 
plain all the symmetry reductions of interest and which 
must be considered in any theory of BFO's structural 
phase transitions: A polar distortion that can in princi- 
ple be oriented along any spatial direction (r^ symme- 
try), and in-phase (M^) and anti-phase (R^) 0$ rota- 
tions around the three Cartesian axes. The atomic dis- 



placements associated with the two AFD order param- 
eters (i.e., the oxygen-octahedra rotations) are uniquely 
defined by symmetry; hence, these modes are trivial in 
this sense. In contrast, the polar distortions are not de- 
termined by symmetry: any combination of -like dis- 
placements of the Bi, Fe, and O sub-lattices is in princi- 
ple valid. Following the usual first-principles approach to 
simple ferroelectric perovskites like BaTi03 or PbTi03,^ 
one would determine the specific atomic displacements 
that define the FE order parameter by computing the 
unstable (soft) polar mode of the cubic phase of the com- 
pound; the result thus obtained for BFO is depicted in 
Fig. 3(a). In materials like BaTi03, such a soft mode 
captures very accurately the atomic distortions associ- 
ated to the relevant FE phases, e.g., tetragonal PAmm 
and rhombohedral Rim. It is not obvious that the same 
will be true for BFO, where we would like to describe 
simultaneously super-tetragonal phases, which imply a 
very large distortion of the cubic cell, and the rhom- 
bohedral ground state, where the polar distortion co- 
exists with very large Oq rotations. Interestingly, we were 
able to demonstrate that the traditional approach works 
well for BFO: We performed a mode-by-mode decompo- 
sition of the atomic distortions connecting the prototype 
PmZm-G phase with the PAmm-C (as representative of 
our large-c/a phases) and i?3c-G structures, and checked 
that the -like component is captured very accurately 
by the soft FE mode of the cubic phase (to within a 93% 
for PAmm-G and 99% for i?3c-G). We can thus conclude 
that it is possible to describe all the FE phases of BFO 
with relatively simple theories that include only one polar 
mode. 

The three primary order parameters described above 
are clearly the driving force for the structural transitions 
in BFO. For a given phase of the material, the occur- 
rence of a particular combination of such primary dis- 
tortions involves a specific breaking of the Pm3m sym- 
metry of the cubic perovskite structure, which in turn 
results in the activation of secondary order parameters 
that become allowed in the low-symmetry phase. The 
most significant secondary distortions that we found in 
our BFO's phases are listed in the last column of Table I 
and sketched in Fig. 3. There is a considerable number 
of such secondary modes; the ones involving the largest 
atomic displacements can be easily grouped in two cate- 
gories: AFE patterns (see (c) to (f) modes in Fig. 3) and 
twisting modes of the 0$ octahedra ((b) in Fig. 3). In 
this sense, BFO is very different from ferroelectrics like 
BaTi03 or PbTi03, where the relevant FE phases do not 
present any secondary modes (note the absence of addi- 
tional distortions for the PAmm, Amm2, and R3m sym- 
metries listed in Table I, which are the relevant ones for 
BaTi03 and PbTiOs). One thus needs to wonder: How 
important are these secondary distortions? Do they play 
a role in determining the energetics and relative stability 
of BFO's phases, or can they be ignored in an effective 
theory of BFO's structural phase transitions';®' 

We quantified the importance of the secondary modes 
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in the following approximate manner: We considered 
the PBEsol+ U equilibrium structures of all the relevant 
phases, artificially set to zero the secondary atomic dis- 
tortions, and computed the energy of the modified struc- 
tures. The obtained energy increments with respect to 
the actual equilibrium phases are very significant: they 
range from tens of meV/f.u. for the monoclinic (Pc-C, 
Cm-C, and Cc-C) phases to more than a hundred for the 
orthorhombic (Pna2\-G and Pnma-G) ones. A more ex- 
act estimate can easily be performed for Pnma-G, as we 
found that in this case the most relevant secondary modes 
are clearly associated to Bi displacements: By fixing the 
Bi ions at their high-symmetry positions and relaxing all 
other structural parameters, we obtained an energy in- 
crease of 125 meV/f.u. with respect to the fully-relaxed 
Pnma-G structure. Thus, our results show that the en- 
ergy changes associated to the secondary modes are of 
the same magnitude as the energy differences between 
different phases, which implies that these modes play a 
key role in determining BFO's phase diagram. In par- 
ticular, the large effects obtained for the orthorhombic 
phases indicate that their stability depends crucially on 
occurrence of the AFE patterns associated to Bi's off- 
centering. We can thus conclude that an effective theory 
of BFO's structural transitions must account for the ef- 
fect of these secondary modes. 



2. Phenomenological theories 

The Devonshire-Landau phenomenological approach 
to phase transitions in bulk ferroelectric s lilliJ anc j its ex- 
tension to epitaxially-constrained films pSHS constitutes 
the simplest, yet powerful, theory that one might try to 
use to model BFO. Working out such a theory for BFO - 
i.e., determining the simplest possible Landau potential 
and temperature dependence of the parameters - consti- 
tutes a great challenge that, as far as we know, remains to 
be tackled. In the following we discuss what our results 
imply as regards the Landau theory of BFO. 

In order to describe all the known phases of this com- 
pound, the corresponding Landau potential should be 
written in terms of a three-dimensional FE polarization 
P (which would correspond to the atomic distortions dis- 
cussed in Section IV.B.l), as well as two AFD order pa- 
rameters associated, respectively, to in-phase and anti- 
phase Og-octahedra rotations. The cross terms between 
these three three-dimensional primary order parameters, 
and the additional terms that will appear if a non-zero 
epitaxial strain is considered, should allow us to repro- 
duce the intricate energy landscape of BFO and its low- 
symmetry minima. 

Indeed, in cases with several order parameters, it is 
possible to obtain stable low-symmetry phases from low- 
order Landau potentials. Imagine, for example, a FE 
perovskite that develops a polarization along the [1,1,1] 
direction as well as an in-phase Og rotation around the 
[0,0,1] axis. Such instabilities can be described with a 




FIG. 6. Energy landscape diagrams as introduced in Ref. 1171 
Filled, open, and shaded symbols correspond to minima, max- 
ima, and saddle points of the energy, respectively. The T, R, 
O, and M labels denote phases with exact tetragonal, rhom- 
bohedral, orthorhombic, and monoclinic symmetries, respec- 
tively, (a) Simplest scenario that gives raise to a monoclinic 
Ma minimum.^ (b) Simplest scenario that gives raise to si- 
multaneous Ma and Mc minima; our discussion focuses on 
the left part of this diagram (see text). 



Landau potential truncated at 4th order in both the FE 
and AFD order parameters. The resulting phase would 
have a monoclinic Pc (Ma) symmetry, exactly as the Pc- 
C structure of Table I. Hence, according to this example, 
it might be possible to describe all BFO's phases with a 
low-order Landau theory. 

However, our results show that the Landau theory for 
BFO would be significantly more complicated, especially 
in what regards the relative stability of the large-c/a 
phases. To illustrate this point, let us consider a sim- 
plified version of BFO in which only FE distortions and 
cell strains are allowed, and try to determine the order 
of the Landau potential F(P) required to describe ferro- 
electricity in such a system. 

In a landmark article, Vanderbilt and CoherP^ ana- 
lyzed the form of the Landau potential needed to describe 
low-symmetry phases in FE perovskites. In essence, they 
showed that a potential F(P) can present tetragonal or 
rhombohedral minima if expanded up to 4th order in P; 
the occurrence of orthorhombic minima requires a 6th- 
order theory, and one needs to go up 8th order to have 
minima of monoclinic symmetry. This work was essential 
to understand which Landau potentials are needed to de- 
scribe the monoclinic phases that were being found at the 
time in perovskite solid solutions such as PbZri^Ti^Oa 
(M A type^and PbZn 1/3 Nb 2/3 3 -PbTi03 [M c type).™ 
The energy landscape associated to an 8th-order po- 
tential with a monoclinic Ma minimum is sketched in 
Fig. 6(a), following the convenient representation scheme 
introduced in Ref. [17] 

We simulated our simplified (FE-only) version of BFO 
by forcing the material to have a 5-atom unit cell in 
which only polar (r^ ) distortions and cell strains are 
allowed. (Of course, this cell was appropriately doubled 
to capture the G- and C-AFM spin arrangements.) If 
we impose such a constraint to the phases in Table I, 
we immediately recover the symmetries that were bro- 
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ken by the AFD modes: The Pc-G and Cc-G phases 
reduce to a single monoclinic Ma structure with space 
group Cm-C; the Cm-C phase changes to a monoclinic 
Mc with Pm-C symmetry; R3c-G gives us a R3m-G 
phase analogous to BaTiCVs ground state, etc. We can 
then consider two additional phases - namely, the super- 
tetragonal PAmm-C and orthorhombic Amm2-G listed 
in Table I -, to sketch the energy landscape of Fig. 6(b). 
(To plot Fig. 6(b), the structural stability against T-like 
distortions of the T and M phases was explicitly checked. 
We have divided the diagram in two sectors to empha- 
size that the distortions connecting the super-tetragonal 
phases with the rhombohcdral and orthorhombic struc- 
tures are very large.) The most notable feature of this 
energy diagram is that it presents two inequivalent mon- 
oclinic minima, as opposed to only one as in Fig. 6(a). 
Further, if we follow the lowest-energy path connecting 
the Ma and Mc minima through triclinic structures, we 
will necessarily cross either a saddle point (case depicted 
in Fig. 6(b)) or a maximum. According to the analysis of 
Ref. H7J the existence of a triclinic saddle point requires 
a Landau potential of 10th order, while a 12th-order the- 
ory is needed to have a triclinic maximum. Note that 
Landau potentials of such a high order are unheard of 
among FE perovskites, even if complex solid solutions are 
considered. Amusingly, in their paper 17 Vanderbilt and 
Cohen justified the interest of discussing theories of very 
high order by writing that "the discovery (or synthesis) 
of a material having such a behavior may be challenging, 
but is by no means impossible." Our analysis shows that 
BFO (even a simplified version of it) is such a material.ESl 



3. Atomistic theories 

Effective theories of the inter-atomic interactions in 
ferroelectric perovskites, with parameters computed from 
first-principles, were introduced in the early 90's by Rabe 
and Vanderbilt. 45 Ever since, these so-called effective 
Hamiltonians have made it possible to perform statistical 
simulations of increasingly complex materials, from crys- 
talline BaTiO^U to disordered PbZri^^Ti^OaPS suc- 
cessfully reproducing temperature-driven phase transi- 
tions, response properties, etc. More recently, an effec- 
tive Hamiltonian for BFO has been derived by Kornev 
et ai.pS who thus extended the approach to incorpo- 
rate magneto-structural interactions in the model. Such 
a groundbreaking development has led to great physi- 
cal insight into BFO's ferroelectric and magnetoelectric 
properties) 26 * 48 ! as well as into the material's behavior 
under applied electrical and magnetic^ fields. On the 
other hand, in view of recent experimental results, we 
now know that some of the model predictions (e.g., the 
occurrence of a 1 4/ 'mem phase at high temperature) are 
questionable. In the following we briefly summarize what 
our results teach us about how to construct an accu- 
rate effective Hamiltonian for BFO, extracting the corre- 
sponding conclusions as regards the theory of Kornev et 



al. 

The first step of the classic approach to constructing 
effective Hamiltonians consists in identifying the relevant 
local distortions that must be retained in the model, so 
that we can use a coarse-grained representation of the 
atoms in the unit cell of our compound. In the case of 
BFO, there are clearly two local distortions that need to 
be considered: (1) a polar displacement pattern compat- 
ible with the FE (I^ ) soft mode of Fig. 3(a), and (2) 
the rotation of individual 0§ octahedra around an ar- 
bitrary axes, whose in-phase (-M3 ) and anti-phase (R~l) 
repetition throughout the crystal reproduces the relevant 
AFD modes. As shown in Section IV.B.l, it is enough to 
consider one local polar mode to reproduce the FE dis- 
tortion of the i?3c-G ground state and large-c/a phases, 
which allows us to work with a relatively simple model. 
A first-principles effective Hamiltonian considering these 
two types of local variables was first constructed to study 
SrTiOs, 23 and this was also the starting point of the work 
of Kornev et al. for BFO. 

In Section IV.B.l we demonstrated the importance of 
the secondary distortions in determining the relative sta- 
bility of BFO's phases. The most relevant secondary 
modes are clearly the Bi-related AFE patterns that oc- 
cur in the Pnma-G and Pna2\-G phases. Fortunately, 
it is possible to incorporate such effects in an effec- 
tive Hamiltonian without extending or complicating the 
model: We can choose the above mentioned local po- 
lar modes to be centered at the Bi atoms, as sketched 
in Fig. 7(a), so that (i) their homogeneous repetition 
throughout the crystal reproduces the FE soft mode of 
Fig. 3(a) and (ii) the zone-boundary modulations repro- 
duce approximately the most relevant AFE distortions of 
the Bi atoms. Note that, alternatively, one could think 
of using local polar modes centered at the Fe atoms (see 
Fig. 7(b)). However, while this option is valid to repro- 
duce BFO's FE distortions, it fails to capture Bi's AFE 
patterns (a zone-boundary modulation of the Fe-centered 
modes results in null Bi displacements) !^ Consequently, 
an effective model based on Fe-centered modes will put a 
considerable energy penalty on the Pnma-G and similar 
phases. Such was the approach adopted by Kornev et 
al, which may explain their prediction that a 1 4 /mem 
structure, and not Pnma, occurs in the phase diagram 
of bulk BFO. 

As regards the rest of (less important) secondary 
modes, it might be possible to incorporate their effect 
by suitably renormalizing the Hamiltonian parameters. 
To make this idea more precise, let us denote by u (resp. 
v) the distortions that will (resp. will not) be explicitly 
considered in the model. The usual effective Hamiltoni- 
ans i? c ff (u), which work well for materials like BaTiOs in 
which secondary distortions are clearly not critical, can 
be formally defined as: 

H eS (u) « E(u,v)\ v=0 , (1) 

where E(u,v) is the first-principles energy of an arbi- 
trary configuration of the compound. In contrast, we 
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(a) (b) 



FIG. 7. Examples of local polar modes that can be used as 
variables of an effective Hamiltonian for BFO. (a) Centered at 
the Bi atom, (b) Centered at the Fe atom. The quantities 81 
identify displacements of the corresponding / atom in the soft 
FE mode of the system; they should be divided by a factor 
that takes into account how many cells share atom /. For 
BFO we obtained S Bi = 0.80, 8 Fc = 0.06, 8 i = -0.42, and 
S 02 = -0.04. 

could define an effective Hamiltonian H c fi(u) designed 
to account for the effect of secondary distortions as: 

H c g(u) min v E(u,v). (2) 

Such a refined approach should improve the accuracy of 
the models in all cases, and it might prove critical to 
obtain correct results for compounds as challenging as 
BFO. The implementation of these ideas remains for fu- 
ture work. 



C. The role of Bismuth 

The Bi cations play a key role in BFO's structural tran- 
sitions. This can be predicted already from very simple 
steric arguments: In BiMC>3 perovskites, where M is a 
first-row transition metal, the lattice parameter is essen- 
tially determined by the ionic radii of the metal and oxy- 
gen ions. This situation, which corresponds to a small 
value of the so-called tolerance factorpS tends to result 
in either the off-centering of the Bi 3+ cation or the occur- 
rence of AFD modes, both of which imply the shortening 
of some Bi-0 bondsPSlThis is exactly what is commonly 
observed in B1MO3 crystals, and the main reason why 
some of these compounds make it possible to combine 
ferroelectricity (related to Bi's off-centering) and mag- 
netism (associated with the transition metals) at high 
temperatures. 

Beyond its relatively small size, Bi 3+ presents an elec- 
tronic configuration (6s 2 p°) that allows for orbital rear- 
rangements suitable to form very directional bonds with 
neighboring oxygen atoms. Such Bi-0 bonds tend to re- 
sult in a lone pair on the non-bonding side, exactly as 
found in BFO's R3c-G phasePThis can be readily visu- 
alized in an electron-localization-function (ELFpSl anal- 
ysis of the computed electronic structure: As shown in 




(c) Pnma-G 



FIG. 8. (Color online.) Electronic-localization-function 
(ELF) maps computed for the _R3c-G, Cc-C, and Pnma-G 
phases. The figures on the left show the isosurface for an 
ELF value of 0.3 superimposed to the atomic structure; on 
the right we show the ELF contour plots in the planes de- 
fined by the labeled ions. We also indicate the shortest Bi-O 
distances (in Angstrom) as obtained from PBEsol+ U calcu- 
lations. 



Fig. 8(a), there is a distinct non-bonding localization do- 
main on the side of Bi that is opposite to the three neigh- 
boring O atoms, which is the signature of a lone pair! 55 * 56 -! 
The occurrence of such a lone pair was discussed at length 
by Ravindran et alW^ on the basis of first-principles cal- 
culations similar to ours; our results for BFO's R3c-G 
phase (Figs. 8(a) and 9(a)) essentially reproduce their 
study.^ 

We computed the ELF maps for the other BFO phases 
found in this work. Figure 8 shows the results for two rep- 
resentative cases: Cc-C and Pnma-G. It is immediate to 
note that a lone pair forms in the super-tetragonal Cc-C 
phase, as might have been expected from Bi's large off- 
centering and the anisotropic spatial distribution of its 
neighboring oxygens. The case of Pnma-G is quite dif- 
ferent, though: As shown in Fig. 8(c), in this phase the Bi 
cations have four neighboring oxygens that form a rather 



12 




Energy (eV) 

(a) _R3c-G 



40,^ — r 




Energy (eV) 

(b) Cc-C 



40 




Energy (eV) 
(c) Pnma-G 

FIG. 9. (Color online.) Electronic density of states for 
the -R3c-G, Cc-C, and Pnma-G phases, as obtained from 
PBEsol+ U calculations. Note that, in these AFM phases, the 
results for the spin-up and spin-down channels are identical. 



regular BiC>4 tetrahedron. The corresponding ELF plots 
show a very isotropic localization domain around Bi. 
There is no clear lone-pair formation in this case; fur- 
ther, such a localization domain is not typical of bond- 
ing electrons, as evidenced by the slightly smaller ELF 
values along the directions of the Bi-0 bonds. Hence, 
it might be more appropriate to interpret this result as 
corresponding to a semi-core-like case. Interestingly, the 
partial density of states results shown in Fig. 9 indicate 
that these three phases are very similar as regards or- 
bital occupation, even if they clearly differ in terms of 
Bi-0 bonding and lone-pair occurrence. Hence, our re- 
sults illustrate Bi's electronic flexibility and its ability to 
form different coordination complexes with the neighbor- 
ing oxygens. 

These chemical effects are clearly the driving force 
for the structural transitions in BFO. Note that all the 
BFO phases discussed here, either ferroelectric or para- 
electric, have an energy that is lower than that of the 
cubic structure by more than 800 meV/f.u. (see Ta- 
ble I). In contrast, the cubic and polar phases differ by 
about 15 meV/f.u. in the case of prototype ferroelectric 



BaTiC>3, where the Coulomb dipole-dipole interactions 
are known to be the driving force for the FE instability.^ 
Noting that BFO and BaTiC>3 are rather similar as re- 
gards the magnitude of the dipole-dipole forces,^ we can 
conclude that such an enormous difference in the strength 
of the structural instabilities must be associated with the 
dominant role of the Bi-0 chemistry in BFO. Then, the 
relative stability of BFO's low-energy phases is probably 
determined by factors that involve smaller energy dif- 
ferences, such as subtle competitions between different 
Bi-0 bonding mechanisms, the build-up of dipole-dipole 
interactions in the FE phases, etc. Analyzing these is- 
sues in detail falls beyond the scope of the present work. 
We hope our findings will stimulate further theoretical 
studies of the chemical bond in these phases, so that 
the factors controlling the occurrence of AFD and/or FE 
distortions (especially the super-tetragonal ones) can be 
elucidated. 

Let us conclude by noting that our results for BFO - 
with most phases being dominated by either AFD or FE 
distortions - are clearly reminiscent of the competition 
between AFD and FE instabilities that is well-known to 
occur in many perovskite oxides. Such a competition has 
been studied in detail in SrTiOapS anc j i s one of the fac- 
tors responsible for the rich phase diagram of materials 
like PbZr 1 _ :c Ti :E 03!sSl Interestingly, BFO is peculiar inas- 
much its FE soft mode in dominated by the A-site cation, 
whereas ferroelectricity in SrTi03 is driven by the B-site 
transition metal and PbZr 1 _ a .Ti :r 03 is an intermediate 
case. Hence, BFO may constitute a new model system 
for the investigation of competing-instability phenomena 
in perovskite oxides. 



V. SUMMARY AND CONCLUSIONS 

We have used first-principles methods to perform a 
systematic search for potentially-stable phases of mul- 
tiferroic BiFc03. We worked with a 40-atom super- 
cell (i.e., a 2x2x2 repetition of the cubic perovskite 
cell) that is compatible with the atomic distortions that 
are most common among transition-metal perovskite ox- 
ides, namely, ferroelectric, anti-ferroelectric, and anti- 
ferrodistortive. We obtained plenty of distinct low-energy 
phases of the compound; here we have restricted the 
discussion to the most stable ones. Many of the ob- 
tained minima present complex structural distortions and 
very low symmetry (e.g., monoclinic Ma and Mc space 
groups) while preserving a relatively small unit cell. As 
far as we know, this is quite unique among perovskite ox- 
ides, as the monoclinic structures reported so far are asso- 
ciated to complex solid solutions (e.g., PbZri-^TLjOs or 
PbZn 1 / 3 Nb2/303-PbTi0 3 ), present large unit cells (e.g., 
BiMn03 and BiScOs), or are obtained under special con- 
ditions (e.g., thin films subject to appropriate epitaxial 
strains or bulk compounds under external electric fields) . 
In contrast, our study shows that bulk BiFe03 presents 
per se a collection of simple low-symmetry minima of the 
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energy. 

Our findings have a number of important implica- 
tions for the research on BiFeC>3 and related materi- 
als. Maybe the most general and interesting one stems 
from the demonstration that BFO can form plenty of 
(meta-)stable structural phases, which suggests that re- 
cent puzzling observations - ranging from possible struc- 
tural transitions at low temperatures^ to surface-s pecifi c 
atomic structures^! and strain-induced new phases^ZI _ 
may just be reflecting BFO's intrinsic structural rich- 
ness. Additionally, our results will provide useful infor- 
mation to the experimental workers exploring the pos- 
sibility of obtaining large functional (piezoelectric, mag- 
netoelectric) effects in BiFe03 films grown on strongly- 
compressive substrates: We have shown that there are 
plenty of phases - all with large polarizations and c/a 
aspect ratios - that can be realized in such conditions, 
including possibilities with monoclinic and orthorhombic 
symmetries. Our results also provide new insights con- 
cerning the relative importance of the various structural 
distortions that can occur in BiFc03, stressing the key 
role that the so-called secondary modes play in determin- 
ing the relative stability of the observed phases. 

Our work also has implications for theoretical stud- 
ies of BiFeOa. First, we present a critical comparison of 
the various DFT schemes most commonly employed to 
study BiFe03 and related compounds, and discuss the 
existing difficulties to quantify the relative phase stabil- 
ity. Second, we draw important conclusions as regards 
the effective modeling of structural phase transitions in 
BiFeOs, in connection with both Landau- type and atom- 
istic theories. Our analysis shows that BiFe03 is rather 
unique, and that its modeling needs to address issues - 
ranging from the work with high-order Landau potentials 



to the accurate treatment of secondary distortions - that 
are unheard of in the work with classic materials such as 
BaTi03, PbZri_2.Ti3.O3, or even relaxor ferroelectrics. 
Finally, our results provide quantitative evidence for the 
dominant role that the Bi-0 bond formation plays in 
BiFeOs's structural instabilities. Further, our analysis 
suggests that some of the phases discussed here do not 
exhibit the "lone-pair mechanism" usually invoked to ex- 
plain the Bi-0 directional bonds in BiFeOs. We take 
this as a new illustration of Bi's ability to form diverse, 
competitive in energy, bonding complexes with its neigh- 
boring oxygens. 

In conclusion, we have used first-principles simulation 
methods to illustrate, quantify, and analyze in some de- 
tail the structural richness of BiFeOs, the most relevant 
representative of the family of Bi-based transition-metal 
perovskite oxides. Our simulations have revealed a vari- 
ety of novel effects, some of which have important impli- 
cations for current experimental and theoretical research 
on this material. We thus hope this work will help clar- 
ify and further stimulate research on these ever surprising 
compounds. 
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